HLA-B*46:01:01:01 and HLA-DRB1*09:01:02:01 are associated with anti-rHuEPO-induced pure red cell aplasia.

Treatment of anemia in patients with chronic kidney disease (CKD) with recombinant human erythropoietin (rHuEPO) can be disrupted by a severe complication, anti-rHuEPO-induced pure red cell aplasia (PRCA). Specific HLA genotypes may have played a role in the high incidence of PRCA in Thai patients (1.7/1,000 patient years vs. 0.03/10,000 patient years in Caucasians). We conducted a case-control study in 157 CKD patients with anti-rHuEPO-induced PRCA and 56 controls. The HLA typing was determined by sequencing using a highly accurate multiplex single-molecule, real-time, long-read sequencing platform. Four analytical models were deployed: Model 1 (additive: accounts for the number of alleles), Model 2 (dominant: accounts for only the presence or absence of alleles), Model 3 (adjusted additive with rHuEPO types) and Model 4 (adjusted dominant with rHuEPO types). HLA-B*46:01:01:01 and DRB1*09:01:02:01 were found to be independent risk markers for anti-rHuEPO-induced PRCA in all models [OR (95%CI), p-values for B*46:01:01:01: 4.58 (1.55-13.51), 0.006; 4.63 (1.56-13.75), 0.006; 5.72 (1.67-19.67), 0.006; and 5.81 (1.68-20.09), 0.005; for DRB1*09:01:02:01: 3.99 (1.28-12.49), 0.017, 4.50 (1.32-15.40), 0.016, 3.42 (1.09-10.74), 0.035, and 3.75 (1.08-13.07), 0.038, in Models 1-4, respectively. HLA-B*46:01:01:01 and DRB1*09:01:02:01 are susceptible alleles for anti-rHuEPO-induced PRCA. These findings support the role of HLA genotyping in helping to monitor patients receiving rHuEPO treatment.

Effects of polymorphisms in the MTHFR gene on 5-FU hematological toxicity and efficacy in Thai colorectal cancer patients.

Background: The two common methylenetetrahydrofolate reductase (MTHFR) polymorphisms 677G>A and 1298A>C may have been affecting 5-FU toxicity in cancer patients for decades. Drug efficacy has also been shown by previous studies to be affected. In this study, we investigated the effects of these polymorphisms on 5-FU hematological toxicity and treatment efficacy, to provide enhanced pharmacological treatment for cancer patients. Methods: This is a retrospective study involving 52 Thai colorectal cancer patients who were treated with 5-FU based therapy, using TaqMAN real-time PCR to genotype the MTHFR polymorphisms (677G>A and 1298A>C). The toxicity and response rate were assessed using standardized measures. Results: Neutropenia was significantly more likely to be experienced (P=0.049, OR=7.286, 95% CI=0.697-76.181) by patients with the MTHFR 677G>A polymorphism, in the same way as leukopenia (P =0.036, OR=3.333, 95%CI=2.183-5.090) and thrombocytopenia (P<0.001, OR=3.917, 95%CI=2.404-6.382). The MTHFR 1298A>C polymorphism had no statistical association with hematological toxicity in 5-FU treatment. The response rate to 5-FU was not significantly affected by these two polymorphisms. Conclusion: The MTHFR polymorphism 677G>A is a significant risk factor for developing leukopenia, neutropenia and thrombocytopenia as toxic effects of 5-FU therapy in cancer patients. Therefore, patients receiving 5-FU-based therapy should be aware of their polymorphisms as one risk factor for experiencing severe toxicity.

Associations between UGT1A1 and SLCO1B1 polymorphisms and susceptibility to neonatal hyperbilirubinemia in Thai population.

Hyperbilirubinemia is the main mechanism that causes neonatal jaundice, and genetics is one of the risk factors of hyperbilirubinemia. Therefore, this study aims to explore the correlation between two genes, UGT1A1 and SLCO1B1, and hyperbilirubinemia in Thai neonates. One hundred thirty seven neonates were recruited from Division of Clinical Chemistry, Ramathibodi Hospital. UGT1A1*28 and *6 were determined by pyrosequencing whereas, SLCO1B1 388A > G and 521 T > C genetic variants were determined by TaqMan® real-time polymerase chain reaction. Neonates carrying with homozygous (AA) and heterozygous (GA) variants in UGT1A1*6 were significantly related to hyperbilirubinemia development compared with wild type (GG; P < 0.001). To the combined of UGT1A1, total bilirubin levels in homozygous variant were higher significantly than heterozygous variant and wild type (P = 0.002, P = 0.003, respectively). Moreover, SLCO1B1 combination was significant differences between the hyperbilirubinemia and the control group (P = 0.041). SLCO1B1 521 T > C variant provide protection for Thai neonatal hyperbilirubinemia (P = 0.041). There are no significant differences in UGT1A1*28 and SLCO1B1 388A > G for the different severity of hyperbilirubinemia. The combined UGT1A1*28 and *6 polymorphism is a strong risk factor for the development of severe hyperbilirubinemia in Thai neonates. Therefore, we suggest neonates with this gene should be closely observed to avoid higher severities of bilirubin.

Genotypic and phenotypic landscapes of 51 pharmacogenes derived from whole-genome sequencing in a Thai population.

Differences in drug responses in individuals are partly due to genetic variations in pharmacogenes, which differ among populations. Here, genome sequencing of 171 unrelated Thai individuals from all regions of Thailand was used to call star alleles of 51 pharmacogenes by Stargazer, determine allele and genotype frequencies, predict phenotype and compare high-impact variant frequencies between Thai and other populations. Three control genes, EGFR, VDR, and RYR1, were used, giving consistent results. Every individual had at least three genes with variant or altered phenotype. Forty of the 51 pharmacogenes had at least one individual with variant or altered phenotype. Moreover, thirteen genes had at least 25% of individuals with variant or altered phenotype including SLCO1B3 (97.08%), CYP3A5 (88.3%), CYP2C19 (60.82%), CYP2A6 (60.2%), SULT1A1 (56.14%), G6PD (54.39%), CYP4B1 (50.00%), CYP2D6 (48.65%), CYP2F1 (46.41%), NAT2 (40.35%), SLCO2B1 (28.95%), UGT1A1 (28.07%), and SLCO1B1 (26.79%). Allele frequencies of high impact variants from our samples were most similar to East Asian. Remarkably, we identified twenty predicted high impact variants which have not previously been reported. Our results provide information that contributes to the implementation of pharmacogenetic testing in Thailand and other Southeast Asian countries, bringing a step closer to personalized medicine.

Host genetic factors of COVID-19 susceptibility and disease severity in a Thai population.

Host genetic factors have been shown to play a role in SARs-CoV-2 infection in diverse populations. However, the genetic landscape differs among various ethnicities; therefore, we explored the host genetic factors associated with COVID-19 disease susceptivity and disease severity in a Thai population. We recruited and genotyped 212 unrelated COVID-19 Thai patients and 36 controls using AxiomTM Human Genotyping SARs-COV-2 array, including 847,384 single nucleotide polymorphisms related to SARs-COV-2 pathogenesis, immune response, and related comorbidity No SNPs passed the genome-wide significance threshold of p value <1 × 10-8. However, with a threshold of p value <1 × 10-5, a locus on chromosome 5q32 was found to have a suggestive association with COVID-19 disease susceptibility (p value 6.9 × 10-6; Q-Q plot λ = 0.805, odds ratio 0.02). Notably, IL17B is a gene located in this linkage disequilibrium block and is previously shown to play a part in inflammation and pneumonia. Additionally, a suggestive locus on chromosome 12q22, harboring EEA1 and LOC643339, was associated with COVID-19 disease severity (p value 1.3 × 10-6 - 4.4 × 10-6, Q-Q plot λ = 0.997, odds ratio 0.28-0.31). EEA1 is involved in viral entry into cells, while LOC643339 is a long non-coding RNA. In summary, our study suggested loci on chromosomes 5q32 and 12q22 to be linked to COVID-19 disease susceptibility and disease severity, respectively. The small sample size of this study may lessen the likelihood that the association found is real, but it could still be true. Further study with a larger cohort is required to confirm these findings.

Influence of SULT1A1*2 Polymorphism on Plasma Efavirenz Concentration in Thai HIV-1 Patients.

PURPOSE: Plasma efavirenz (EFV) concentrations within therapeutic levels are essential to successfully treat patients suffering from human immunodeficiency virus (HIV) type 1. In addition to the drug-metabolizing enzyme CYP2B6, other phase II drug-metabolizing enzymes and transporters may have an important role in the pharmacokinetics of EFV. Thus, the influence of phase II drug-metabolizing enzymes and drug transporters on plasma EFV levels was investigated in Thai HIV patients receiving EFV. PATIENTS AND METHODS: Genotyping was performed by TaqMan® real-time PCR in 149 HIV-infected Thai adults, and plasma efavirenz concentration was measured by a validated high-performance liquid chromatography in 12 hours after dosing steady-state plasma samples at week 12 and 24. RESULTS: Patients with three or more copies of SULT1A1 had significantly lower median plasma EFV concentrations than those carrying two copies at week 12 (p=0.046) and SULT1A1*2 (c.638G>A) carriers had significantly lower median plasma EFV concentrations compared to those not carrying the variant at week 24 (p=0.048). However, no significant association was found after adjusting for CYP2B6 genotype. CONCLUSION: Genetic variation in a combination of SULT1A1*2 and SULT1A1 copy number may contribute to variability in EFV metabolism and thereby may impact drug response. The influence of a combination between the SULT1A1 and CYP2B6 genotype on EFV pharmacokinetics should be further investigated in a larger study population.

Phenotype prediction and characterization of 25 pharmacogenes in Thais from whole genome sequencing for clinical implementation.

Publicly available pharmacogenomics (PGx) databases enable translation of genotype data into clinically actionable information. As variation within pharmacogenes is population-specific, this study investigated the spectrum of 25 clinically relevant pharmacogenes in the Thai population (n = 291) from whole genome sequencing. The bioinformatics tool Stargazer was used for phenotype prediction, through assignment of alleles and detection of structural variation. Known and unreported potentially deleterious PGx variants were identified. Over 25% of Thais carried a high-risk diplotype in CYP3A5, CYP2C19, CYP2D6, NAT2, SLCO1B1, and UGT1A1. CYP2D6 structural variants accounted for 83.8% of all high-risk diplotypes. Of 39 known PGx variants identified, six variants associated with adverse drug reactions were common. Allele frequencies of CYP3A5*3 (rs776746), CYP2B6*6 (rs2279343), and NAT2 (rs1041983) were significantly higher in Thais than East-Asian and global populations. 121 unreported variants had potential to exert clinical impact, majority were rare and population-specific, with 60.3% of variants absent from gnomAD database. This study demonstrates the population-specific variation in clinically relevant pharmacogenes, the importance of CYP2D6 structural variation detection in the Thai population, and potential of unreported variants in explaining drug response. These findings are essential in development of dosing guidelines, PGx testing, clinical trials, and drugs.

CYP2D6 Predicts Plasma Donepezil Concentrations in a Cohort of Thai Patients with Mild to Moderate Dementia.

PURPOSE: Donepezil, a drug frequently used to treat dementia, is mainly metabolized by cytochrome P450 2D6 (CYP2D6). This study investigated the relationships between CYP2D6 genotype and activity scores as well as predicted phenotype of plasma donepezil concentrations in 86 Thai dementia participants. MATERIALS AND METHODS: CYP2D6 was genotyped using bead-chip technology (Luminex xTAG® v.3). Steady-state trough plasma donepezil concentrations were measured using high-performance liquid chromatography. RESULTS: Sixteen genotypes were found but the most frequent genotypes detected among our participants were CYP2D6*10/*10 (27.9%) and *1/*10 (26.7%). One-third of the participants had an activity score of 1.25 which predicted that they were normal metabolizers. The overall median (interquartile range) of plasma donepezil concentration was 51.20 (32.59-87.24) ng/mL. Normal metabolizers (NMs) had lower plasma donepezil concentrations compared to intermediate metabolizers (IMs) (41.15 (28.44-67.65) ng/mL vs 61.95 (35.25-97.00) ng/mL). Multivariate analysis showed that CYP2D6 activity score (r2 = 0.50) and the predicted phenotype (independent of dose) could predict the plasma donepezil concentration (r2 = 0.49). CONCLUSION: Plasma donepezil concentration in NMs was lower compared to IMs. Additional studies with larger sample size and use of next-generation sequencing as well as its outcomes are warranted to confirm the benefit of using pharmacogenetic-guided treatment for donepezil.

SLCO1B1 and ABCG2 Gene Polymorphisms in a Thai Population.

INTRODUCTION: Genetic polymorphisms of drug transporters influence drug transporter activity and alter pharmacokinetic profiles of the drugs. Organic anion transporting polypeptide 1B1 (OATP1B1) and breast cancer resistance protein (BCRP) are important transporters encoded by solute carrier organic anion transporter family member 1B1 (SLCO1B1) gene and ATP-binding cassette subfamily G member 2 (ABCG2) gene, respectively. Polymorphisms in these genes are associated with increased plasma statins concentrations, statin-induced myopathy and poor response to allopurinol treatment. PURPOSE: We explored allele and genotype frequencies of SLCO1B1 and ABCG2 genes including their predicted phenotypes in 53 Thai participants. Of these, 17 had chronic kidney disease and were on statins. MATERIALS AND METHODS: Genotyping analysis for SLCO1B1 c.521T>C (rs4149056), c.388A>G (rs2306283), g.-11187G>A (rs4149015), and ABCG2 c.421C>A (rs2231142) was done by using TaqMan® Real time PCR. All were tested for Hardy-Weinberg Equilibrium. RESULTS: Most of the participants (80%) had normal function haplotypes SLCO1B1 (*1A and *1B) while decreased (*5, *15, and *17) and unknown (*21) function haplotypes were less observed. Four phenotypes of SLCO1B1 were observed: 69.81% had normal function (*1A/*1A,*1A/*1B, and *1B/*1B), 13.21% had intermediate function (*1A/*17, *1B/*15 and *1B/*17), 9.43% had indeterminate function (*1A/*21 and *1B/*21) and 7.55% had low function (*5/*15, *15/*15, and *15/*17). ABCG2 c.421A allele frequency was 25%. The frequency of ABCG2 c.421CA and AA phenotypes were 37.7% and 5.7%, respectively. The allele and genotype frequencies observed are consistent with reports in Asians. However, there were differences in major allele distributions between Asians and Caucasians for SLCO1B1 c.388A>G; SLCO1B1 c.388G were highly found in Asians, but c.388A were more observed in Caucasians. CONCLUSION: This study showed that in the Thai population, there were 4 SNPs of SLCO1B1 and ABCG2 genes. This finding may be clinically applied to minimize inter-individual variability of drugs such as statins and allopurinol. Further study with a larger sample size is needed to assess the drug profiles and responses to treatment.

Effect of drug metabolizing enzymes and transporters in Thai colorectal cancer patients treated with irinotecan-based chemotherapy.

Genetic polymorphisms in drug metabolizing enzymes and drug transporters may affect irinotecan toxicity. Although genetic polymorphisms have been shown to influence the irinotecan toxicity, data are limited in Thai population. Thus, the aim of this study was to assess the allele and genotype frequencies and the relationship between CYP3A4/5, DPYD, UGT1A1, ABCB1, and ABCC2 genetic variations and irinotecan-induced toxicity in Thai colorectal cancer patients. One hundred and thirty-two patients were genotyped, and the effect of genetic variations on irinotecan-induced toxicity was assessed in 66 patients who received irinotecan-based chemotherapy. Allele frequencies of ABCB1 c.1236C > T, ABCB1 c.3435C > T, ABCC2 c.3972C > T, ABCG2 c.421C > A, CYP3A4*1B, CYP3A4*18, CYP3A5*3, DPYD*5, UGT1A1*28, and UGT1A1*6 were 0.67, 0.43, 0.23, 0.27, 0.01, 0.02, 0.64, 0.19, 0.16, and 0.09, respectively. DPYD*2A and DPYD c.1774C > T variants were not detected in our study population. The ABCC2 c.3972C > T was significantly associated with grade 1-4 neutropenia (P < 0.012) at the first cycle. Patients carrying both UGT1A1*28 and *6 were significantly associated with severe neutropenia at the first (P < 0.001) and second (P = 0.017) cycles. In addition, patients carrying UG1A1*28 and *6 had significantly lower absolute neutrophil count (ANC) nadir at first (P < 0.001) and second (P = 0.001) cycles. This finding suggests that UGT1A1*28, *6, and ABCC2 c.3972C > T might be an important predictor for irinotecan-induced severe neutropenia.

Pharmacogene Variation in Thai Plasmodium vivax Relapse Patients Treated with a Combination of Primaquine and Chloroquine.

PURPOSE: Pharmacogenes have an influence on biotransformation pathway and clinical outcome of primaquine and chloroquine which are often prescribed to treat Plasmodium vivax infection. Genetic variation may impact enzyme activity and/or transporter function and thereby contribute to relapse. The aim of the study was to assess allele, genotype frequencies and the association between pharmacogenes variation and primaquine response in Thai patients infected with Plasmodium vivax. PATIENTS AND METHODS: Fifty-one patients were genotyped for 74 variants in 18 genes by Sequenom MassARRAY® and Taqman® SNP Real-Time PCR. RESULTS: SNP frequencies were not significantly different between relapse (n=4) and non-relapse (n=47) patients. However, the CYP2C19 c.681G>A, the frequency of the A-allele that defines the non-functional CYP2C19*2 haplotype was significantly higher compared to the G-allele (OR=5.14, p=0.021). Patients heterozygous for ABCG2 c.421C>A had a higher odds ratio (OR=8.75, p=0.071) and the frequency of the G-allele of UGT2B7 c.372G>A was higher compared to the A-allele (OR=3.75, p=0.081). CYP2C19, ABCG2 and UGT2B7 emerged as potential high priority genes. CONCLUSION: Decreased activity of CYP2C19, ABCG2 and UGT2B7 in combination with CYP2D6 intermediate or poor metabolizer status may expose patients to a higher risk of Plasmodium vivax relapse. Further investigations are warranted to substantiate these findings.

CYP2D6 genotype analysis of a Thai population: platform comparison.

The highly polymorphic CYP2D6 gene locus leads to a wide range of enzyme activity. Since there are limited data for Thai, the major aim was to investigate CYP2D6 genetic variation in a large Thai population (n = 920). CYP2D6 genotyping was performed using four different platforms. Genotype call rates of the Luminex xTAG® and AmpliChip CYP450 test were 96.5% and 87.4%, respectively. Based on Luminex xTAG® data, the most common alleles and genotypes were *1 0 (49.6%), *1 (24.6%), *2 (10.8%), *5 (6.7%), *41 (6.5%) and *1/*10 (23.9%), *10/*10 (21.5%), *2/*10 (9.4%), *5/*10 (6.9%), *10/*41 (5.7%), respectively. Challenges and limitations of the platforms evaluated are discussed. These data add to our knowledge regarding interethnic variability in CYP2D6 activity and contribute to improving drug therapy in the Thai.

UGT1A1 polymorphisms associated with prolactin response in risperidone-treated children and adolescents with autism spectrum disorder.

The aim of this study was to investigate the association of drug-metabolizing enzyme and transporter (DMET) polymorphisms with the risperidone-induced prolactin response using an overlapping gene model between serum prolactin level and hyperprolactinemia in autism spectrum disorder (ASD) patients. Eighty-four ASD patients who were receiving risperidone for at least 1 month were recruited and then assigned to either the normal prolactin group or the hyperprolactinemia group based on their serum prolactin level. The genotype profile of 1936 (1931 single nucleotide polymorphisms (SNPs) and 5 copy number variation (CNVs) drug metabolism markers was obtained using the Affymetrix DMET Plus GeneChip microarray platform. Genotypes of SNPs used to test the accuracy of DMET genotype profiling were determined using TaqMan SNP Genotyping Assay kits. Eighty-four patients were selected for the allelic association study after microarray analyses (51 in the normal prolactin group, and 33 in the hyperprolactinemia group). An overlapping allelic association analysis of both analyses discovered five DMET SNPs with a suggestive association (P < 0.05) with risperidone-induced prolactin response. Three UGT1A1 SNPs (UGT1A1*80c.-364C > T, UGT1A1*93 c.-3156G > A, and UGT1A1 c.-2950A > G, showed a suggestive association with the risperidone-induced prolactin response and found to be in complete linkage disequilibrium (D' value of 1). In this DMET microarray platform, we found three UGT1A1 variants with suggestive evidences of association with the risperidone-induced prolactin response both measured by hyperprolactinemia and by prolactin level. However, due to the lack of validation studies confirmation and further exploration are needed in future pharmacogenomic studies.

Pharmacogenomics and Efficacy of Risperidone Long-Term Treatment in Thai Autistic Children and Adolescents.

The purpose of this study was to evaluate the association of pharmacogenomic factors and clinical outcome in autistic children and adolescents who were treated with risperidone for long periods. Eighty-two autistic subjects diagnosed with DSM-IV and who were treated with risperidone for more than 1 year were recruited. Pharmacogenomics and clinical outcome (CGI-I, aggressive, overactivity and repetitive score) were evaluated. Almost all patients showed stable symptoms on aggressive behaviour (89.02%), overactivity (71.95%), repetitive (70.89%) behaviour and all clinical symptoms (81.71%). Only 4.48% of patients showed minimally worse CGI-I score. Patients in the non-stable symptom group had DRD2 Taq1A non-wild-type (TT and CT) frequencies higher than the clinically stable group (p = 0.04), whereas other gene polymorphisms showed no significant association. Haplotype ACCTCAT (rs6311, rs1045642, rs1128503, rs1800497, rs4436578, rs1799978, rs6280) showed a significant association with non-stable clinical outcome (χ2  = 6.642, p = 0.010). Risperidone levels showed no association with any clinical outcome. On the other hand, risperidone dose, 9-OH risperidone levels and prolactin levels were significantly higher in the non-stable compared to the stable symptom group (p = 0.013, p = 0.044, p = 0.030). Increased appetite was the most common adverse drug reaction and associated with higher body-weight, whereas it was not significantly associated with genetic variations and non-genetic information. In conclusion, risperidone showed efficacy to control autism, especially aggressive symptoms in long-term treatment. However, Taq1A T - carrier of dopamine 2 receptor gene - is associated with non-stable response in risperidone-treated patients. This study supports pharmacogenomics testing for personalized therapy with risperidone in autistic children and adolescents.

CYP2D6 polymorphisms and their influence on risperidone treatment.

Cytochrome P450 enzyme especially CYP2D6 plays a major role in biotransformation. The interindividual variations of treatment response and toxicity are influenced by the polymorphisms of this enzyme. This review emphasizes the effect of CYP2D6 polymorphisms in risperidone treatment in terms of basic knowledge, pharmacogenetics, effectiveness, adverse events, and clinical practice. Although the previous studies showed different results, the effective responses in risperidone treatment depend on the CYP2D6 polymorphisms. Several studies suggested that CYP2D6 polymorphisms were associated with plasma concentration of risperidone, 9-hydroxyrisperidone, and active moiety but did not impact on clinical outcomes. In addition, CYP2D6 poor metabolizer showed more serious adverse events such as weight gain and prolactin than other predicted phenotype groups. The knowledge of pharmacogenomics of CYP2D6 in risperidone treatment is increasing, and it can be used for the development of personalized medication in term of genetic-based dose recommendation. Moreover, the effects of many factors in risperidone treatment are still being investigated. Both the CYP2D6 genotyping and therapeutic drug monitoring are the important steps to complement the genetic-based risperidone treatment.

Pharmacogenomic Study Reveals New Variants of Drug Metabolizing Enzyme and Transporter Genes Associated with Steady-State Plasma Concentrations of Risperidone and 9-Hydroxyrisperidone in Thai Autism Spectrum Disorder Patients.

The present study sought to investigate the genetic variants in drug metabolizing enzyme and transporter (DMET) genes associated with steady-state plasma concentrations of risperidone among Thai autism spectrum disorder (ASD) patients. ASD patients taking risperidone for at least 1 month were enrolled for this pharmacogenomic study. Genotyping profile was obtained using Affymetrix DMET Plus array interrogating 1931 variants in 231 genes. Steady-state plasma risperidone and 9-hydroxyrisperidone were measured using liquid chromatography/tandem mass spectrometry assay. The final analysis included 483 markers for 167 genes. Six variants, ABCB11 (c.3084A > G, c.∗420A > G, c.∗368G > A, and c.∗236G > A) and ADH7 (c.690G > A and c.-5360G > A), were found to be associated with plasma concentrations of risperidone. 9-Hydroxyrisperidone and the total active-moiety levels were associated with six gene variants, SCLO1B1 (c.-11187G > A and c.521T > C), SLCO1B3 (c.334G > T, c.699A > G, and c.1557G > A), and SLC7A5 c.∗438C > G. Polymorphisms in UGT2B4 c.∗448A > G and CYP2D6 (c.1661G > C, c.4180G > C, and c.-2178G > A) showed considerable but not significant associations with metabolic ratio. This pharmacogenomic study identifies new genetic variants of DMET genes in monitoring risperidone therapy.